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We discuss the non-linear extension to the integrated Sachs-Wolfe effect (ISW) resulting from 
the divergence of the large scale structure momentum density field. The non-linear ISW effect leads 
to an increase in the total ISW contribution by roughly two orders of magnitude at I ~ 1000. This 
increase, however, is still below the cosmic variance limit of the primary anisotropies; at further small 
angular scales, secondary effects such as gravitational lensing and the kinetic Sunyaev-Zel'dovich 
(SZ) effect dominates the non-linear ISW power spectrum. We show this second-order non-linear 
ISW contribution is effectively same as the contribution previously described as a lensing effect due 
to the transverse motion of gravitational lenses and well known as the Kaiser-Stebbins effect under 
the context of cosmic strings. Due to geometrical considerations, there is no significant three point 
correlation function, or a bispectrum, between the linear ISW effect and its non-linear extension. 
Q_C The non-linear ISW contribution can be potentially used as a probe of the transverse velocity of dark 

matter halos such as galaxy clusters. Due to the small contribution to temperature fluctuations, 
of order few tenths of micro Kelvin, however, extracting useful measurements on velocities will be 
challenging. 
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I. INTRODUCTION 

> 

\q | It is by now well known the importance of cosmic microwave background (CMB) temperature fluctuations as a 
probe of cosmology Qj. The accuracy to which cosmological information can be extracted depends on how well we 
understand individual processes that lead to anisotropies in CMB temperature. Though effects during recombination 
are now well understood , contributions and modifications to CMB anisotropies due to large scale structure between 
last scattering surface and today is not completely well established. This is primarily due to the non-linear evolution of 
the large scale structure at low redshifts such that simple analytical calculations based on linear theory may no longer 
& ■ be applicable. In general, large scale structure affects CMB through two processes: gravity and Compton scattering. 
Or The modifications due to gravity arises from frequency changes via gravitational red and blue-shifts, while during the 
reionized epoch, photos can both generate and erase primary fluctuations through scattering via free electrons. 

Here, we discuss an effect due to gravitational redshift commonly known in the literature as the integrated Sachs- 
Wolfe (ISW; ||]) effect at late times.Q The temperature fluctuations in the ISW effect results from the differential 
redshift effect from photons climbing in and out of time evolving potential perturbations from last scattering surface to 
^ , present day. In currently popular cold dark matter cosmologies with a cosmological constant, significant contributions 
J*"! ■ arises at redshifts less than 1 and on and above the scale of the horizon at the time of decay. 

rS | Here, in particular, we extend previous discussions on the ISW effect, usually due to linear fluctuations in the 
density field, to the non-linear regime of clustering. The non-linear contribution to ISW effect is generally called 
the Rees-Sciama (RS; Q) effect, though such a distinction is arbitrary since temperature fluctuations in CMB arise 
essentially due to time evolving potentials in both the linear and non-linear regime of fluctuations. Using the continuity 
equation, we show that the non-linear contributions to ISW effect comes from the divergence of the large scale structure 
momentum density field. We model the large scale structure momentum density field using the recently popular halo 
model where large scale structure density field fluctuations can be described through dark matter within halos and 
correlations between halos m. Our analytical calculations are consistent with those of jjj based on numerical 
simulations. 

In addition to a recalculation of the contribution, we also discuss this non-linear extension to the ISW effect under 
the context of proposed contributions to CMB temperature in the literature. We show the non-linear ISW contribution 
is essentially same as temperature anisotropy produced through the transverse motions of foreground gravitational 
lenses as was first discussed in Ref. M and with further discussions and corrections in Ref . M . The effect is generally 



To avoid confusion, we distinguish contributions to ISW effect during matter dominant era as late ISW effect, though there 
is an additional contribution during the radiation dominance. 
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called the moving lens contribution [[10 1 1 and its cosmic string analogue is the well known Kaiser-Stebbins effect |L2| . 
Under this lensing description, the effect can be described as the gravitational lensing of the CMB dipole, formed by 
the transverse motion of a halo or a galaxy cluster, in the rest frame the halo. We show that there is no significant 
non-Gaussian correlations associated with the non-linear ISW effect and the bispectrum formed with linear effects are 
negligible. There is also no correlation between the kinetic SZ effect, due to the momentum density field along the 
line of sight, and the non-linear ISW effect. 

The layout of the paper is as follows. In § y, we outline the non- linear contribution to the ISW effect after reviewing 
briefly background material relevant for current calculations in the context of the adiabatic cold dark matter (CDM) 



models. In § III, we discuss our results and consider the non-linear contribution under the context of suggested 
secondary effects in CMB. We also discuss the bispectrum formed by the combined linear ISW and its non-linear 
extension. In the same section, we study the possibility for an extraction of the transverse velocity of galaxy clusters 
using the non-linear ISW effect. 

II. CALCULATIONAL METHOD 

We first review the properties of adiabatic CDM models relevant to the present calculations. 

A. Adiabatic CDM Model 

The expansion rate for adiabatic CDM cosmological models with a cosmological constant is 

H 2 = Hi [n m (i + zf + Q K (1 + zf + n A ] , (1) 

where Hq can be written as the inverse Hubble distance today Hq 1 = 2997.9/i _1 Mpc. We follow the conventions 
that in units of the critical density 3Hq/8ttG, the contribution of each component is denoted f^, i = c for the CDM, 
g for the baryons, A for the cosmological constant. We also define the auxiliary quantities £l m — £l c + fif, and 
= 1 — f2j, which represent the matter density and the contribution of spatial curvature to the expansion rate 
respectively. 

Convenient measures of distance and time include the conformal distance (or lookback time) from the observer at 
redshift z = 

r dz' 

and the analogous angular diameter distance 

d A = fl^fi" 1 / 2 sinh(£T ^ /2 r) . (3) 

Note that as £Ik — 0, cLa — > r and we define r[z = oo) = tq. 

In linear theory, the density field may be scaled backwards to higher redshift by the use of the growth function 
GO), where S(k,r) = G(r)5(k,0) g| 

G(r)oc^) r dz'{l + z')(^\ . (4) 



Ho J z(r) v '\H(z') 

Note that in the matter dominated epoch G tx a = (l + z) -1 . 

Although we maintain generality in all derivations, we illustrate our results with the currently favored ACDM 
cosmological model. The parameters for this model are ft c = 0.30, ilf, = 0.05, £!a = 0.65, h — 0.65, Y p — 0.24, 
n = 1, X = 1, with a normalization such that mass fluctuations on the 8h Mpc -1 scale is as = 0.9, consistent 
with observations on the abundance of galaxy clusters |L4j and COBE normalization |L5|. A reasonable value here is 
important since higher order correlations is nonlinearly dependent on the amplitude of the density field. To compute 
the linear power spectrum, we adopt the fitting formula for the transfer function given in [061 . 



2 



B. ISW effect 



The integrated Sachs- Wolfe effect results from the late time decay of gravitational potential fluctuations. The 
resulting temperature fluctuations in the CMB can be written as 



T ISW (n) 



dr<b(r, nr) . 



(5) 



where the overdot represent the derivative with respect to conformal distance (or equivalently look-back time). Writing 
multipole moments of the temperature fluctuation field T(n), 

(6) 



<„,„ = J dnT(n)>T*(n), 
we can formulate the angular power spectrum as 

For the ISW effect, multipole moments are 

TSW _ A I a K 



(7) 



— 2 I dr$(k)/ ; (fc)^ m (k), 



(8) 



with Ii{k) = J drW lsw (k,r)ji(kr), and the window function for the ISW effect, W lsw (see below). The angular 
power spectrum is then given by 



where the three-dimensional power spectrum of the time-evolving potential fluctuations are defined as 

<$(kl)$(k 2 )> = (27T) 3 fo(k 1 +k 2 )Pii(fcl). 



(9) 



(10) 



The above expression for the angular power spectrum can be evaluated efficiently under the Limber approximation 
,7j for sufficiently high I values, usually in the order of few tens, as 



r isw _ / . \w_ J_ p . . 

°i — I ar d 2 



k = — ,r 
d A 



(11) 



In order to calculate the power spectrum of time-derivative of potential fluctuations, we make use of the cosmological 
Poisson equation |18|. In Fourier space, we can relate the fluctuations in the potential to the density field as: 



(12) 



Thus, the derivative of the potential can be related to a derivative of the density field and the scale factor a. Considering 
a flat universe with SIk — 0, we can write the full expression for the power spectrum of time-evolving potential 
fluctuations, as necessary for the ISW effect valid in all regimes of density fluctuations, as 



9 / £l m \ 2 ( Hn x 



P*i(*,r)= t — ^ 



P M (k,r)-2^P st (k,r)+P st (k,r) 



(13) 



with VF ISW = -2 in equations | and 0. 

To calculate the power spectrum involving the correlations between time derivatives of density fluctuations, Py, 
and the cross-correlation term involving the density and time-derivative of the density fields, P S g , we make use of the 
continuity equation, which in position, or real, space can be written in the form: 
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5(x,r) = -V- [l + <J(x,r)]v(x,r). (14) 

In the linear regime of fluctuations, when <5(x, r) = G(r)J(x, 0) <C 1, the time derivative is simply 5 lln (x, r) = 
— V • v(x, r) and we can obtain the well-known result for linear theory velocity fluctuations in Fourier space as 



iG8{k,Q)-^. (15) 



k 

'k 

Thus, in linear theory, P^ = k 2 P vv (k,r) = G 2 P]' 5 n (fc,0) and P gS = kP Sv {k,r) = GGP^(k,0) 
These lead to the well-known results for the linear ISW effect, with a power spectrum for 4> as 



lin/; \ _ ^ / fl m \ ( Hq 



** V ' ' 4 V a J V k 



n 2 



--G(r) + G 



a 



Plf(k,0). (16) 



The term within the square bracket is F 2 where F = G/a following derivation for the linear ISW effect in |19| . Even 
though, we have replaced the divergence of the velocity field with a time-derivative of the growth function, it should 
be understood that the contributions to the ISW effect comes from the divergence of the velocity field and not directly 
from the density field. Thus, to some extent, even the linear ISW effect reflects statistical properties of the large scale 
structure velocities. 

In the mildly non- linear to fully non- linear regime of fluctuations, the approximation in equation [l4|, involving 
5 <C 1, is no longer valid and a full calculation of the time-derivative of density perturbations is required. This can 
be achieved in the second order perturbation theory, though, such an approximation need not be fully applicable 
as the second order perturbation theory fails to describe even the weakly non-linear regime of fluctuations exactly. 
Motivated by applications of the halo approach to large scale structure and results from numerical simulations 
p?H , we consider a description for the time-derivative of density fluctuations and rewrite equation |l4| as 

<j(x, r) = - V ■ v(x, r) - V • <5(x, r)v(x, r) , (17) 

where we have separated the momentum term involving p = (1 + S)v to a velocity contribution and a density velocity 
product. In Fourier space, 



5(k) =ik-p(k) 

d 3 k 



^k.v(k)+ J — 5(k-k')ik.v(k'), 



(18) 



where we have dropped the time-dependence for clarity. The first term involving the velocity field leads to the linear 
theory ISW effect, while the non-linear aspects are captured in the term involving convolution of the Sv term (see, also 
0). Following the approach motivated by Hu in Ref. discussed by Cooray in Ref. |(| and investigated in detail 
through numerical simulations by Sheth et al. in Ref. |20[| , we can write the power spectrum of density derivatives, 
or equivalently the divergence of the momentum density field, as 

■fjj(^) = k 2 P pp {k) 

( r7 3 k' 



k2 S W 1 ^^'^''^ 

d 3 k' f d 3 k" 

rH'fi"Tssw(k. - k', -k - k", k', k" 



(2tt) 3 J (2tt) 3 ' 

(19) 

In the non-linear regime, k k' ~ k and k — k" ~ k and we can simplify by integrating over angles to obtain, 

i r k' 2 rlh' 

P„ = k 2 P^{k)+'-k 2 P ss (k) J 1^-P vv (k'). 

(20) 

Note that in the deeply non-linear regime, contribution from the trispectrum formed by the velocity-density corre- 
lations drop out; this is due to the fact that under the halo approach, the non- linear trispectrum resulting from the 
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single halo term is independent of the configuration and thus, not on // and (see, discussions in |6j and |22|). Also, 
the term involving product of cross-power spectra between the density and velocity fields does not contribute in the 
non-linear regime; at small scales, the density fluctuations are independent of the large scale velocity field. In figure p], 
we show the power spectrum of momentum divergence, which in linear theory is simply described by the divergence 
of the velocity field, with the extension to the non-linear regime of fluctuations following equation 
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Our description of the non-linear momentum power spectrum is similar to the derivation of the non-linear kinetic 
SZ effect in Hu |Ilj and Cooray || (see, also, |2(|) which involves the momentum field along the line of sight. In 
equation the integral over the velocity power spectrum is simply the RMS of the velocity fluctuations 

<n S = J dk^f^ . (21) 

Thus, the non-linear power spectrum of the momentum field divergence, as relevant for the ISW effect, involves one 
of three components of the velocity field with 1 /3w^ ms ; note that the expression for the fully non-linear contribution 
to the momentum field along the line of sight is similar to above and also involves one component of the velocity 
field as discussed in pl[ | and |(|. The resulting expression for the non-linear momentum field is fully consistent with 
simulations pfj. 

In addition to the power spectrum of density derivatives, in equation |l3|, we also require the cross power spectrum 
between density derivatives and density field itself P s $. In M), using the halo approach as a description of the 
momentum density field observed in simulations, it was found that the cross-correlation between the density field and 
the momentum field can be well described as 

P P s (k) = y/ ' P pp (k)P 5S (k) . (22) 

This is equivalent to the statement that the density and momentum density fields are perfectly correlated with a 
cross-correlation coefficient of 1; this relation is exact at mildly-linear scales while at deeply non-linear scales this 
perfect cross-correlation requires mass independent peculiar velocity for individual halos [ pot . Using this observation, 
we make the assumption that P S g ~ -JPssP^, which is generally reproduced under the halo model description of 
the cross-correlation between density field and density field derivatives. This cross-term leads to a 10% reduction of 
power at multipoles between 100 and 1000, when compared to the total when linear and non-linear contributions are 
simply added. 



III. DISCUSSION 

In figure ||, we show the angular power spectrum of the ISW effect with its non-linear extension (which we have 
labeled RS for Rees-Sciama effect). The curve labeled ISW effect is the simple linear theory calculation with a power 
spectrum for potential derivatives given in equation |16i The curves labeled "lin" and "nl" shows the full non-linear 
calculation following the description given in equation]l3| and using the linear theory or full non-linear power spectrum, 
in equation [2^, for the density field, respectively. For the non-linear density field power spectrum, we make use of 
the halo approach for large scale structure clustering |^|,^| and calculate the power spectrum through a distribution 
of dark matter halos. We use linear theory to describe the velocity field in both linear and non-linear cases; since the 
velocity field only contributes as an overall normalization, through t> rms , its non-linear effects, usually at high k values, 
are not important due to the shape of the velocity power spectrum and the behavior of the integral in equation p]l 

As shown in figure the overall correction due to the non-linear ISW effect leads roughly two orders of magnitude 
increase in power at I ~ 1000. The difference between linear and non-linear theory density field power spectrum in 
equation only leads to at most an order of magnitude change in power. Note that the curve labeled "lin" agrees 
with previous second order perturbation theory calculations of the Rees-Sciama effect M, while the curve labeled "nl" 
is also consistent with previous estimates based on results from numerical simulations. 

There is an additional feature that should be observed in figure |J but not properly described previously in the 
literature. At I ~ 100 to 1000, there is roughly a 10% decrease in total power from what is generally described in the 
literature as the Rees-Sciama contribution. This is due to the P$ s term in equation [l3|. If one simply adds the linear 
ISW and non-linear RS contributions, this cross-correlation term is not present. This dip, due to a cancellation, is 
present when comparing results based on perturbation theory for RS and numerical simulations for the full non-linear 
ISW effect (see, 0); The cross term provides a natural explanation for the slight decrease in power. 

In figure @, we break the contribution to the non-linear effect, without the linear ISW or cross-term contributions, 
as a function of steps in redshift. As shown, one essentially finds equal contributions over a wide range in redshift 
with most of the contributions coming from a redshift ~ 1 to multipoles of few hundred where the power peaks. 
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Even though there is roughly two orders of magnitude increase in power at multipoles around 1000, in figure |4|, we 
show that this increase is still below the cosmic variance associated with primary contribution given by 

ao =\/Sf c < cmb - (23) 

where A; is the bin size in multipole space and / s k y is the fraction of sky covered. The signal-to-noise for the detection 
of the power spectrum is 




where C" is the power spectrum of noise with C™ = Cp MB + Cf + Cf et ; Cf is any contribution from secondary effects 
and Cf et is any detector noise contribution. As shown in figure |J(b), the cumulative signal-to-noise is less than one 
even for a full sky experiment with / s k y = 1 and no detector noise. Here, we include contributions from the thermal 
and kinetic Sunyaev-Zel'dovich effects (SZ; |2^||) and lensing effect on CMB [||J as secondary contributions to the 
noise power spectrum. The top line with cumulative signal-to-noise slightly above one is when thermal SZ effect is 
removed from the noise contribution. 

Given the significant sample variance and the fact that the non-linear ISW effect has no special property, such as 
a different frequency spectrum from thermal CMB as in the case of thermal SZ effect |2^,^], it is unlikely that the 
non-linear ISW effect power spectrum can be extracted from CMB data easily. Going to much smaller angular scales, 
or large multipoles, leads to a reduction in the cosmic variance, though, increase in power associated with other small 
angular scale secondary effects such as the thermal SZ effect or the kinetic SZ effect can complicate any detection of 
the non- linear ISW effect. Later, we will address the possibility whether, instead of statistical properties such as the 
power spectrum or the bispectrum (see below), we can extract the associated signal from individual objects such as 
massive galaxy clusters. 



A. ISW Bispectrum 

Following our earlier discussion for the angular power spectrum, we can also consider the bispectrum (see, [ p6[ ) , or 
the Fourier analog of the three point correlation function: 

B(n,m,i) = (T(fi)T(m)T(i)) (25) 

where the sum is over (li, TOi), (fa, TII2), (h-, 1713). Statistical isotropy again allows us to express the correlation in terms 
of an rn-independent function, 

(ai imi ai 2 m 2 ai 3 7n 3 ) = ( ^ ^ B hhh . (26) 

Here the quantity in parentheses is the Wigner-3j symbol. The orthonormality relation for Wigner-3j symbol implies 

h h h 

m\m 2 rn 3 

For the coupling between ISW effects, using multipolar moments written in equation ^, we can write the bispectrum 

as 



B hhh= X! ( m \ m 2 m 3 ) ( a hr ni ai 2m2 ai 3m3 ) . (27) 



(2Ii + l)(2I 2 + l)(2/ 3 + l) h h h 



Bllhh V 4tt ^000 

x b hihil3 , (28) 

with 



G 



°h,h,l 3 



— J k\dk\ J k\dk 2 J k^dk 3 



x / x 2 dxj h (k 1 x)ji 2 (k 2 x)ji 3 (k 3 x) . 

As before, Limber approximation |l7| ] allows one to simplify the integrals for speedy calculation 

W lsw (h , r) W lsw (k 2 , r) W lsw (k 3 , r) 



(29) 



dr- 
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x B iii (k 1 ,k 2 ,k 3 ) ]jr , , 



(30) 



Similar to the power spectrum in equation 10, the three dimensional bispectrum of the derivatives of potential 
fluctuations are defined as 



($(fc 1 )$(fc 2 )$ nl (fc 3 )> - (2 7 r) 3 < 5 D (k 1 + k 2 + k 3 ) 
x B^(k x ,k 2 ,k 3 ) 



(31) 



As an approximation, this three-dimensional bispectrum can be calculated in second order perturbation theory (e.g., 
p7[). To obtain an exact result valid both in the weakly non-linear and non- linear regimes, we consider the coupling 
between two linear ISW effects and the non-linear extension involving the V • 5v term. This leads to a bispectrum 
following: 



B. 



x GG 



P ss (k u 0)Ps v (k 2 ,0) 



27 fn n 

a 

ks-k 



H 6 



h 2 h 2 h 2 



Perm. 



-G + G 



(32) 



where permutations are with respect to the ordering of ki,k2 and k 3 leading to a total of six terms. In equation pQ, 
W ISW = —2 as defined earlier. Because of the dependence on an angle, say k 3 • k 2 in B^^(ki,k 2 ,k 3 ), there is 
significant cancellations and the final projected angular bispectrum of the ISW effect is smaller than a simple order of 
magnitude estimate involving the cube of the temperature fluctuation amplitude. In figure [7], we show the cumulative 
signal-to-noise for the ISW bispectrum, with 



— /sky 



R2 

E n hhh 

h,h,h ll h 13 



(33) 



where again the noise power spectrum is C™ (see, Ref. |n| for details). In figure ^, we take C™ = Cf MB to consider 
sample variance from primary anisotropics alone. The cumulative signal-to-noise is in the order of few times 10 -7 , 
suggesting that the bispectrum is unlikely to be detected; similar values are also found for other three point statistics 
such as the skewness or the third moment. Thus, consistent with the second order perturbation theory result [p7| , 
there is no significant non-Gaussian signal at the three-point level formed by the correlation between linear ISW 
effects and the non-linear ISW effect in the non-linear regime of fluctuations. Also note that there is no bispectrum 
of the form (<b(ki)<t nl (k 2 )<b nl (k 3 )) since such a term leads to an odd number of velocity or density fluctuation terms. 



B. Cross-Correlations 



An additional way to extract or detect the presence of this signal is through cross-correlations with another source 
of anisotropy or a tracer of large scale structure. Since the momentum density field is involved in the non-linear ISW 
effect, one can expect the presence of a correlation between another tracer of the momentum density field. It is well 
known that the kinetic SZ effect traces the line of sight large scale structure momentum density field such that the 
temperature fluctuations can be written as a modulation of the velocity field by baryon fluctuations 

T kSZ (n) = / drg{r)n ■ v(r, hr)6 b (r, nr) , (34) 
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where St, is the fluctuation in the baryon field and g(r) is a weight function for Compton-scattering with g = re T = 
X7J 0.0691(1 - Y p )Q g h(l + z) 2 e~ T where r(r) = ^ drf is the optical depth out to r 
The correlation between the non-linear ISW effect and the kinetic SZ effect is then 

(35) 

and is equal to zero through the angular terms involving the integral over k • k' = fc/c'/z'; the cross-correlation involves 
a cosine term from the divergence of the momentum associated with non-linear ISW effect and a sine term from the 
line of sight momentum associated with the kinetic SZ effect. The geometrical cancellation is merely a statement that 
though locally line of sight velocity field may be correlated with its transverse component, on average over the whole 
sky, there is no such correlation. This cancellation can be avoided by several techniques. In we discussed a similar 
situation involving the cross-correlation between the thermal SZ effect and the kinetic SZ effect and suggested the 
use of a quadratic correlation involving the square of the density field. This is also equivalent to the use of absolute 
values of the temperature fluctuations. 



C. Relation to other effects 

We now suggest the non-linear extension to the ISW effect is essentially the same contribution described by || 
using transverse motion of lcnsing objects; under the context of cosmic strings, this contribution is well known as the 
Kaiser-Stebbins effect Q. Writing, T(n) = / drAT(rh, r), we consider the ISW effect such that AT(rn, r) = -2$(r). 
Following our earlier discussion related to the non-linear contribution, we can relate the potential fluctuations to the 
density fluctuations using the Poisson equation. The contributions now follow as AT oc [— a/aS + 5}. The first term 
containing the time-derivative of the scale factor was recently reintroduced by [M as a time-delay effect. As discussed 
in |50|, the time-delay contribution to CMB is second order as it involves a product of the spatial gradient of CMB 
at the last scattering surface and the cumulative time-delay contribution and is not simply described by a first order 
contribution. Additionally, as discussed for the linear ISW effect, only considering the a term leads to an overestimate 
of the fluctuation as there is a cancellation from a first order term involving the S term. 

As before, using the complete continuity equation, we can relate the time derivative of the density fluctuations, 
S, to the divergence of the velocity field and the divergence of the product of density overdensity and velocity 
AT oc [a/aS + V • (1 + 5)v] Here, we consider the non-linear contribution resulting from the AT nl oc V • <5v term and 
reintroducing the density field in terms of potentials using the Poisson equation, 

AT nl = -2V • ($v) 
w -vl • (2V r $) 

= — v sina<5icn cos <fi ■ (36) 

The simplifications assume that potential fluctuations are embedded in a velocity field with much larger coherence 
scale so that gradients in the velocity field do not contribute to temperature anisotropics. Furthermore, we have 
introduced the lensing deflection angle <5i cn = 2V r $ where the gradient is now an angular gradient on the sky and 
there is no contribution to temperature anisotropy from the gradient of potential along the line of sight. This forces 
the contributing component of the velocity field to be the one on the sky and not along the line of sight. Here, a is 
the angle between the line of sight and the velocity field and <f> 1S the position angle from the observer. As it is now 
clear, this latter description of the non-linear ISW effect is what has been provided elsewhere under the context of 
moving gravitational lenses . 

The correspondence between the two effects can also be noted using the the description which is well utilized to 
calculate the gravitational lensing effect on CMB |2j,|l^,^l|. Writing, 

T(n) = T(n + An) 

= T p (n)+V r T p (n)-An, (37) 

where angular deflections due to gravitational lensing follows such that An = 2V r </>, where the projected lensing 
potential is 0(n) = Di s /D s §. Here, T p (n) is the primary CMB contribution while V r T p (n) is the angular gradient on 
the sky of these temperature fluctuations. Essentially, gravitational lensing angular deflections remap the distribution 
of temperature fluctuations and this remapping is captured by the dependence on distances where distance from lens 
to source (CMB last scattering surface) is Di s and the distance from observer to source is D s . Also, as written, 
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gravitational lensing effect on CMB is second order due to the dependence on the angular gradient of the CMB on 
the sky, V r T; this is consistent with the fact that lensing does not change the surface brightness and only results in 
a modification of the temperature fluctuation distribution. 

The second order effect related to non-linear 1SW effect is due to the lensing of the dipole created by the motion 
of halos such that V r T = wsina. To be consistent with the description given under the time-derivative of the 
density fluctuations, we require Di s — D s , such that there is no dependence on the ratio of distance factors, usually 
encountered in lensing studies. This is equivalent to the case that Di <C D s such that Di s ~ D s , which is not 
necessarily true even for clusters at z ~ 0.5. In either case, the non-linear ISW effect cannot be considered as a 
remapping of the temperature fluctuations similar to the conventional gravitational lensing, as there is no real source, 
a temperature gradient, in the problem; the temperature gradient essentially exists from the coordinate transformation 



from a moving lens to a stationary lens 
when attempting to describe the contri 
lensing effect. 



^ The inconsistency with lensing mapping description arises, unfortunately, 
Dution from the divergence of the momentum density field as a gravitational 



D. Towards transverse velocities 



As discussed before, the detection of the power spectrum of temperature anisotropies due to the non-linear ISW effect 
is heavily affected by the dominant cosmic variance of CMB primary anisotropies. Due to geometrical considerations, 
there is also no significant non-Gaussian contribution or a cross-correlation with other effects. Thus, instead of 
statistical properties such as the power spectrum, one can try to extract the associated signal in individual objects 
such as massive galaxy clusters. 

In figure 0, we show the distinct signature formed by the non-linear ISW effect for a cluster of mass 5 x 10 14 Mq 
with a transverse velocity of 100 km sec -1 across the line of sight. We assume a NFW profile for the dark matter 
distribution of the cluster and take a description for the concentration mass relation following fl33| . In terms of the 
density distribution of the cluster, the deflection angle, at an impact distance of r\ from the cluster center, is given by 

h*M = ^ /" r±dr± Y>x) (38) 
c V Jo 



^(r x ) = / dr || p(r±,r \\ ). (39) 



where the surface mass density is 



In above, r± is the distance across the line of sight, r« is the distance along the line of sight and r v is the virial radius 
of the cluster, which we take to be at an overdensity of 200 following the NFW ^] description. Above equation for 
the deflection angle also assumes a circularly-symmetric dark matter distribution within the halo. 

As shown in figure ||, the contribution to temperature fluctuations are at most 0.3 /iK. In order to detect this small 
signal, one has to be able to extract it from other contributions to CMB due to galaxy clusters. The most significant 
contributions from clusters arise from the SZ thermal effect ^3) due to inverse-Compton scattering of photons via 
hot electrons. Here, temperature changes of order 1 mK is produced and these are now routinely observed towards 
massive clusters p3|. The SZ thermal effect, however, has a distinct frequency dependence and in multifrequency 
CMB data, the effect can be separated out from thermal CMB and other fluctuations [E5[|. The next significant 
contribution comes from kinetic SZ effect due to the line of sight motion of the scatterers in clusters (^3|. In figure 0, 
we show the kinetic SZ effect for the same cluster as in figure |5j. Here, we have assumed that the electron distribution 
in clusters is described by the hydrostatic equilibrium || and have taken a line of sight velocity of 100 km sec -1 . The 
contribution due to the SZ kinetic effect, and also SZ thermal effect, is highly peaked towards the center of the cluster 
and can be as high as few tens /xK. 

The next important contribution is due to gravitational lensing of the large scale CMB gradient, whose rms is of 
order 13 fiK arcmin -1 . With a deflection angle of order ~ 0.5 arcmins, the contribution due to lensing is in the 
range of few fiK. As shown in figure || the lensing effect has the same profile shape as the contribution resulting 
from the transverse velocity; the two profiles need not lie in the same direction since the large scale CMB gradient 



2 We note that this issue has led to some confusion in the literature when calculating the so-called moving lens effect with 
some including the ratio of Di s /D s 0] while others, correctly, not (H]] 
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and the transverse velocities can be aligned differently. When the thermal SZ effect is separated, in temperature 
fluctuations, a galaxy cluster exhibits a slight offseted dipolar pattern with a significant temperature increment or 
decrement towards the center, resulting from the direction of the line of sight velocity associated with the kinetic SZ 
effect. Detecting such a profile will certainly remain a challenging goal for the future cluster observations |p5| . 

Eventually, if the transverse velocity contribution can be detected, its amplitude and the direction of the dipole pat- 
tern, provides significant information on large scale structure velocities not generally available from other observations. 
We can ask how well one can detect a typical cluster through this effect. First, to make a reasonable signal-to-noise 
detection, it is clear that one must extract effects such as the kinetic SZ and CMB lensing. This requires detailed 
knowledge on the baryon and dark matter distributions within clusters. Assuming such information is available, we 
can obtain an estimate for the signal-to-noise by noting that the observed signal can be written as |35|]36| l 

^L(0) = S (9)+n(8), (40) 

where s(9) is the profile of the signal and n{9) is the profile of the noise distribution. In order to remove the excess 
noise associated with large scale temperature fluctuations from primary anisotropies or other secondary effects, we 
can construct an appropriately normalized filter which provides an optimal detection of the signal in the presence of 
such noise. This filter can be written as 



= ^ 



d 2 l \s(l)\- 
{2-kY C? 



(41) 



where C™ is the power spectrum of noise. The signal to noise for the detection of the profile, for an axi-symmetric 
distribution for s(6), is 

sy_ fidu{if 



N/ J 2vr C™ 

In figure g we show the cumulative signal-to-noise, as a function of I. Here, we assume a noise power spectrum given 
by the sum of the intrinsic CMB and secondary effects, where secondary effects include all thermal contributions: 
CMB lensing, kinetic SZ and a model for the inhomogeneous reionization. At small is', corresponding to outer extent 
of the cluster, fluctuations in the intrinsic CMB temperature confuse the detection of the signal while at large I's, 
corresponding to inner extent of the cluster, fluctuations from local universe complicate. Even if detailed properties 
of clusters are known so that intrinsic CMB lensing and the kinetic SZ effect can be perfectly separated from the 
transverse effect, it is unlikely that we will not know all sources of temperature fluctuations along the line of sight 
towards a given cluster leading to a confusion in the detection. For typical transverse velocities of order few 100 
km sec -1 , we find signal-to-noise values of order 0.1 suggesting that detection of this signal for individual clusters 
will remain challenging even for an experiment with no instrumental noise contributions; alternatively, one can put 
an upper limit on the transverse velocity contribution at the level of few thousand km sec -1 ; this upper limit is 
considerably larger than what one expects under currently popular ACDM cosmological models. 

The signal-to- noise for detection may be improved if cross-correlation techniques can be considered, for e.g., if one 
has some knowledge on the intrinsic and secondary fluctuations towards the observed cluster and has some knowledge 
on the direction of the large scale bulk flows from other methods, even if amplitude of that bulk flow is not known. 
Separately, since the coherence scale of bulk flows are much larger than an individual cluster, it may be possible to 
extract the transverse effect by averaging the signal over a number of clusters. Using numerically simulated cluster 
images, and realistic sources of noise and confusions, we plan to study how well such an extraction can be performed 
in a future paper. 



IV. SUMMARY 



We have discussed the non-linear extension to the integrated Sachs- Wolfe effect (ISW) resulting from the divergence 
of the large scale structure momentum density field. This non-linear ISW effect, calculated under the recently popular 
halo approach to non-linear large scale structure clustering, leads to an increase in the total ISW contribution by 
roughly two orders of magnitude at I ~ 1000. This increase, however, is still below the cosmic variance limit of the 
primary anisotropies; at further small angular scales, secondary effects such as gravitational lensing and the kinetic 
Sunyaev-Zel'dovich (SZ) effect dominates the non-linear ISW power spectrum. 

Further, we have shown that this second-order non-linear ISW contribution is effectively same as the contribution 
previously described as a lensing effect due to the transverse motion of gravitational lenses and well known as the 
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Kaiser-Stebbins effect related to cosmic strings. Due to geometrical considerations, there is no significant three point 
correlation function, or a bispectrum, between the linear ISW effect and its non-linear counterpart. The correlation 
between the non-linear ISW effect and the kinetic SZ effect is again zero due to geometry associated with the line of 
sight and divergence of the momentum density field. The non-linear ISW contribution can be potentially used as a 
probe of the transverse velocity of dark matter halos such as galaxy clusters, however, due to the small contribution to 
temperature fluctuations of order few tenths of [iK, extracting useful measurements on velocities will be challenging. 
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contribution while the curve labeled "lin" is the contribution resulting from the momentum field under the second order 
perturbation theory. 
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FIG. 3. The redshift dependence of the non-linear ISW contribution. Here, we break contributions at z less than 0.5 (dotted 
line), between 0.5 and 1.0 (dashed line), between 1.0 and 2.0 (long dashed line), and for redshifts greater than 2 (dot-dashed 
line). As shown, most of the contributions arise at a redshift of ~ 1, leading to the peak in power at / ~ 100 and a sharp 
decrease thereafter. 





FIG. 4. (Top): The total ISW contribution compared with cosmic variance due to primary anisotropies. As shown, con- 
tributions are far below the cosmic variance and the non-linearities are not likely to impact the upcoming measurements of 
cosmological parameters from CMB data. At small angular scales, where the non-linear contributions are well above the 
sample variance, other secondary effects, such as the kinetic SZ, dominate the temperature fluctuations, (bottom) cumulative 
signal-to-noise for the detection of the non-linear ISW power spectrum, in the presence of the SZ thermal effect (dotted line) 
and not (solid line), and for an all sky experiment with no detector noise contribution. 
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FIG. 5. The contribution due to the non-linear ISW effect for a cluster of mass 5 x 10 Mq with a transverse velocity of 
100 km sec -1 across the line of sight. The temperature fluctuations produce a distinct dipolar pattern on the sky and are of 
order few tenths of micro Kelvin. Here, x and y coordinates are in terms of the scale radius of the cluster, based on a NFW 
profile. 
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FIG. 6. Temperature fluctuations due to clusters. Top: kinetic SZ effect, middle: lensing of CMB primary temperature 
fluctuations, and bottom: the total contribution from kinetic SZ, lensing and transverse velocities. We use the same cluster as 
shown in figure H. 
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FIG. 7. The cumulative signal-to-noise for the detection of the non-linear ISW effect bispectrum as a function of the multipole 
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FIG. 8. The signal-to-noise for the detection of the non-linear ISW effect. The galaxy cluster considered here is the one in 
figure B 
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